#include<iostream>
using namespace std;
#include<fstream>
#include<cmath>

#include"library.h"

double u(double alpha, double r){
	return r*exp(-alpha*r);
}

double EL(double alpha, double r){
	return -0.5*alpha*alpha + (alpha-1)/r;
}

int main(){

	ofstream fp("vmcHydrogen_aveEL.txt");
	for(double alpha=0.8; alpha <=1.2; alpha += 0.01){
		fp << alpha << "\t" << vmc1d_aveEL(alpha, 100, u, EL, 0, 5, 0.2, 10000) <<endl;
	}
	fp.close();
	return 0;
}
